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Abstract: We compute the virtual 0{a^ nj) corrections to the heavy quark vector current 
correlator in terms of expansions in the external momentum and as an exact numerical 
solution. As a byproduct, the available high-energy expansion at the three-loop level is 
extended. 
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1. Introduction 

Correlators of two currents, basic objects in Quantum Field Theory, provide important 
information for both theoretical and phenomenological applications. Depending on the 
Lorentz structure of the current under consideration, different interesting observables are 
directly related to these quantities, like the hadronic cross section in electron-positron 
annihilation, R{s), and decay rates of Z- and Higgs-bosons. 

On this account current correlators are investigated thoroughly in perturbation theory, 
where even high order calculations are possible. Some applications, like the determination 
of the charm and bottom quark masses through sum-rules carried out to four loops in |2[, 
corrections to the p parameter |^ ^ or the recent evaluation of Os from low energy data 
at the same level of precision |5| , require asymptotic expansions of heavy quark correlators 
in the low and/or high energy limits. Other studies, however, necessitate the knowledge 
of the full external momentum dependence p'^, for instance the determination of the fine 
structure constant at the Z-boson mass scale, aem{Mz) 

Up to the three-loop level, all physically relevant correlators have been computed including 
the full quark mass dependence in 0, |8|, ^ by deriving Fade approximants from asymp- 
totic expansions. At four-loop accuracy, in the case of the vector current including one 
heavy quark, the first two terms in the low energy series are currently known [l^ and 
were obtained by direct Taylor expansion of all propagator-type integrals and subsequent 
reduction of tadpole diagrams. In the high energy limit, besides the leading massless con- 
tribution |n|, 0, , the first two mass correction terms of the absorptive part 



of the scalar, vector and axial- vector current correlators are also available [0, |l^] , partially 



even to five loops [|l9|, . 

Unfortunately, the techniques used to achieve these results are, due to their huge computa- 
tional complexity, not suitable for the computation of higher order terms in the expansions. 
Thus, reconstructing the full momentum dependence at 0{a^) in an analogous manner to 



- 1 - 



the three-loop case demands a completely different method. A promising and effective 



approach in this context is based on differential equations, originally proposed in |21, 22 1. 
Using this technique, we were already able to compute the first 30 terms in the low energy 
expanded polarization function at 0{a1) in |23]. Recently, these results were confirmed 



in utilizing the same approach. 

The purpose of this paper is to exploit this method and demonstrate its applicability up 
to the four-loop level. We compute the vacuum polarization for low and high energies 
as expansions up to {p'^)^^ and respectively, and moreover the full momentum 

dependence in numerical form in the Euclidean and Minkowskian region. As new results, 
we extend the available information in the high energy domain at 0{a^) and provide the 
double fermionic contributions to the vector current correlator at (^(af). 
This work is organized as follows. In the next section we sum up the relevant notation. 
Calculational methods are briefly presented in Section 3, whereas the first terms of the 
expansions and the numerical results are given in Section 4. Finally, we give our conclusions. 

2. Definitions 

Given the vector current j'^{x) = qhix)7^Qhix) composed of the heavy quark field qh{x) 
with mass m, the two point correlator is defined by 

U^'^p^)=i J dxe'''P{0\Tf{x)f{0)\0) (2.1) 

where p'^ is the external momentum. A convenient representation of the tensor 11'*'' (p^) 
through the scalar vacuum polarization function n(p^) is given by 

n'*-(/) = (^_p2gt^u ^ pMp-)n(/) + pt^p-ULip^). (2.2) 

Transversality of the vector current correlator requires Yliip^) = 0. Through the optical 
theorem the aforementioned hadronic ratio R{s) is related to the current correlator 

R(s) = f ^ hadrons) ^ ;^27rlmn(/ = s + ie). (2.3) 
a{e+e~ Ai+Ai'j 

In the framework of perturbation theory, the polarization function can be expanded in the 
strong coupling as 

n(i'=) = <3^ll?E(^)'n"='(P^). (2.4) 

fc>0 ^ ^ 

Since in this work the focus is set on double fermionic contributions at 0{a^) stemming 
form diagrams sketched in fig. |l], it is convenient to decompose n(^)(p^) into bosonic and 
fermionic contributions, 



n(3)(p2) = c'^n^l\p^) + c'^CAU^^i,ip') + CFClu^^l,(j^') (2.5) 
+c2 Tp ng^,(p2) + Cf Ca Tf ngLb') 
+CpT|n(?(p2), 
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Figure 1: Diagrams contributing to n^£ (p^). Solid (dashed) lines refer to massive (massless) 
propagators. 



where Ilj,^ (p^) and Xl^f (p^) denote terms proportional to nj and nj, respectively, with nj 
being the number of active flavours. Cp refers to the Casimir operator and Tp to the trace 
of the fundamental representation of SU{N). Because of the mass hierarchy in the quark 
sector only the heavy quark is considered massive whereas all lighter (n/ = n/ — 1) quarks 
are treated as massless. 

Using the fact that the double fermionic contribution contains no massless cuts, H^^^ (p"^) 
can be expanded in a simple power series in the low energy limit — > 0. In the large 
energy limit p'^ — > —oo, however, non-integer powers of p'^ arise and lead to additional 
logarithms of the form log(— p'^/m(^)^). Thus, defining z = p^ / Am{ii)^ , we end up with 
the expansions 

n>0 

5^C-^(M)^-"log-(-z). (2.7) 

n,m 



3. Calculation 

The basic idea for computing the full p^-dependence of the vacuum polarization function 
is to deal with massive propagator-type integrals. Although the number of integrals is 
moderate (approx. 10^ for the whole four-loop contribution) they pose a challenge as 
far as the reduction to a small set of master integrals is concerned, since two variables, 
z and d = 4 — 2e, the dimension of space-time, are involved. For the determination of 
these master integrals, however, an efficient method through differential equations exists 
and allows for asymptotic expansions to high orders in the external momentum and high 
precision numerics. 

In a first step all Feynman diagrams contributing to the double fermionic corrections have 
been projected onto scalar integrals, which were subsequently reduced to a set of 46 mas- 



ter integrals with the help of integration- by-parts (IBP) identities |25| and the Laporta 



algorithm pq] implemented in IdSolver [2?]. 



To obtain the master integrals, we used the scaling property of propagator-type integrals 
Pi{p^,m^) 
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to get the characteristic differential equation 

.2 d 



(3.1) 



/^P*(/,m2) = -nl^-^p^{p^^rr^)+D[P,{p^,rr^)]Pi{p^,rr^) (3.2) 



with D[Pi{p'^, rn?)] being the mass dimension of Pi{p'^,m'^). Using relations generated from 
the reduction, the right hand side can again be expressed through master integrals which 
leads to a coupled system of inhomogeneous differential equations 



^P^iz) 
az 



Aij{z,e)Pj{z). 



(3.3) 



The matrix Aij{z, e) is composed of rational functions of z and e. Its block-triangular form 
simplifies the problem to a set of several small coupled systems of differential equations 
and therefore provides a systematic approach to solutions of Pi{z) in terms of expansions 
in the external momentum z (or y = 



In the low- (z — > 0) and high-energy (y 0) limits, the system was solved by ansatze similar 
to eqns. ( p.6| , |2.7D for each master integral and coefficients of the series were determined 
recursively up to high powers in z and y, respectively. 

Boundary conditions in the low energy limit are given by massive tadpole diagrams depicted 
in fig. |^(a), analytically calculated in [28|. In the opposite limit, all boundary conditions 
were determined from automatized diagrammatic large momentum expansions which lead 
to products of at most three loop massive tadpoles with massless propagators. Fig. §(b) 
illustrates the needed propagators at the four-loop level. 

The approximate linear complexity of this procedure allows for the computation of coeffi- 
cients to, at least in principle, arbitrary depths. In this work, we concentrate on the first 
30 and 15 terms of the expansions for z < 1 and y < 1, respectively, and are thus able to 
compute the master integrals in those regions with high precision. 




(b) 



Figure 2: Boundary conditions appearing in the expansions of master integrals 



In addition, an exact numerical solution of eqn. (|3.3| ) was accomplished. For this purpose, 
the e-expanded system of differential equations was directly integrated by means of the 
FORTRAN package ODEPACK |29], using the high precision values at a starting point 
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Figure 3: Integration contour chosen for the numerical integration 

l^l ^ 1. Due to the absence of thresholds in the Euchdean domain, the numerical integra- 
tion can easily be carried out along the real axis for z < 0. In the Minkowskian half-plane, 
however, (pseudo)thresholds occur and integration along the real axis is only possible below 
these special points. On the other hand, by virtue of contour deformation into the complex 
plane (see fig. |^, the master integrals can be solved numerically for arbitrary values of 
z > 0, even above z = 1. Variation of the integration contour is used to estimate the real 
achieved precision. Furthermore, the Mellin-Barnes (MB) method has been used at a few 
points to check the obtained values from direct integration. For this purpose, MB repre- 
sentations have been automatically generated with the package MBrepresentation pO|| , 
analytically continued in e and numerically integrated with help of the MB package 
We observed that the high-energy expansion works very well above at least z > 5 and 
therefore the numerical integration was undertaken up to 2 = 10. This fact is also used as 
a cross-check of the high-energy expansion against the numerically computed polarization 
function. 

In a last step, we performed the renormalization of the mass m{fi), the strong coupling 
as(/i) and the external current in the MS-scheme. 

4. Results 

All analytical results for the coefficients in the small- and high-energy series in the MS- 
scheme up to the 30th and 15th term, respectively, are available in Mathematica format 
together with the source of this paper from http : / / arxiv . org. Apart from the new 
expressions at the four-loop level, we also provide so far unknown terms needed for renor- 
malization in the large-energy expansion at the three-loop level. The Taylor series for 
z — > is already known to sufficient depth p3| , p^ ]. 

For the sake of clarity, here we give only the first five terms of both series at 0{a^). The 
first coefficients in the limit z — > are given by 

_ 163868 3287 203 14483 236 . 8 3 / 262877 

^ " "98415" " 2430 + 324 '"^^ " 21870 ™ ^ V262440 
116 29 ,4 1421 4 38909 203 3779 , 

486 ^2 + ^^2vr +^^^ " 19440 ^'+ 324 ^"^^'" 21870 
472 ,2 i^rsA 2/^ 42173 112^ 1784 ^36_ 2 8 ,3 ^ 
+ 3645 ^ 135 '"J + V32805 " 135 ^' + 3645 + 3645 '"+135 '"^ 

(4.1) 
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_n 1842464707 2744471 , 14203 , , 676663 , 1468 16 ,o 

/^U /■ _|_ 7 /- / lA _|_ 7o 

^ ~ 646652160 1064448 27648 870912 42525 315 

/ 95040709 2029 2029 ,0 2 2029 ,4 99421 4 12159109 , 

+ n; 04 H ^nvr 1% ^ TT Cs 

' V 62705664 1728 41472 ^ 41472 ^ 4976640 4644864 ^'^ 

14203 ^ _ 153560977 ^ _ 2936 ^2 32 ^3 \ ^2 / 15441973 _ 32 

"^27648 ''^ ~ 326592000 ™ ~ 42525 315 "7 V 19136250 ~ 45 ^ 

195679 _ 1468 2 16 3 \ 

^637875 ~ 42525 ^ 315 ^ ^ ^ 

_ 56877138427 _ 6184964549 12355 _ 224445289 _ 310736 2 
^ ~ 12609717120 ~ 1556755200 20736 ~ 244944000 ~ 4465125 

128 ,3 / 60361465477 1765 1765 ,4 1765 . o 17297 4 

H C + 04 H Zovr'' H TT^ 

2835 V 29393280000 1296 31104 ^ 31104 ^ 746496 

57669161 12355 37155186557 621472 2 256 3 \ 
— . - Cs + x^r^:r:r tm Cs — , , Im. , . ^- t»r, + t»r. 



iZ5 ,0 / DU0Di4D04<' ^ i/DO i/DO ,4 i^DO ,9 9 i/^ 

C + "-i a'4 ^9 H ^971" H 

2835 V 29393280000 1296 31104 ^ 31104 ^ 746 

57669161 , 12355 , , 37155186557 , 621472 ,o 256 ,3 

A _|_ ] t / / — u / 

17418240 ^ 20736 ""^^ 60011280000 "" 4465125 2835 
2^ 31556642272 256 15480824 310736 2 128 3 \ 

^ V 49228003125 405 ''^ 52093125 4465125 ^ 2835 ^ ' ' 

_ 270605350139987 _ 692437613459 2522821 4400766606253 

^ ~ 40351094784000 ~ 119558799360 ^ 3538944 ~ 4138573824000 

1082216 ^2 256 3 ^ / 432564184014463 360403 ^ 360403 ^4 

~ 12629925 6237 + V 165542952960000 ~ 221184 ~ 5308416 ^ 

360403 ,2 2 17659747 4 44387709491 , 2522821 , ^ 143778598477 , 
5308416 ^ 637009920 10899947520 3538944 ""^^ 189621927936 

_ 2164432 ^2 512 3 \ 667234795424 512 

" 12629925 " ^ 6237 ""y ^ V 1253204308125 ~ 891 ''^ 
1213878224 _ 1082216 2 ^56_ 3 \ 

"^3978426375 ~ 12629925 "" 6237 ""y 

1089264797862809 13766684598973 , 1239683 , , 

/^U A _|_ 7 /- 

^ ~ 114328101888000 1693749657600 ^^1474560 ""^^ 



1089264797862809 _ 13766684598973 1239683 _ 9879141778777 

114328101888000 1693749657600 1474560 ''^ ~ 8136640512000 
1064635136 2 1024 3 ^ / 2680576384798219 1239683 1239683 4 

"11345882625 27027 ™ + V 838813667328000 ~ 645120 ""^ ~ 15482880 ^ 

1239683 ,2 2 8677781 . 1942577613953 1239683 , , 

H H tt' C3 H Cs 

15482880 ^ 265420800 398529331200 ^ 1474560 ^ 

_ 580030100117579693 _ 2129270272 ^2 2048 ^3 \ 
"644112110315520000 ~ 11345882625 ™ 27027 ""J 
^2 / 61443753281463008 _ 2048 160295587080064 ^ _ 1064635136 ^2 

V 136221219619340625 ~ 3861 511075282843125 ~ 11345882625 
1024 .\ , , 

where = log(m(^)^//i^), I2 = log(2), = Lij(l/2) with Lij(x) the polylogarithm 
function and Q the Riemann zeta numbers. 
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^ 5832 324 ^'^ 4320 9 ^ 3 36 ^ 36 ^ 24 

324 "^"^ 27 ^ 27 ^ 216 108 9 ^ 27 27 ' 

_2 2 , 2 2 A 2 22327 358 ^ 20 ^ 19 



224 27 27 216 108 9 ' 

_2 2 , 2 2 A 2 22327 358 ^ 20 

g 'sm 'i;/* + 27 9 =^7* / '^^ \ 2916 81 9 

_J_ 2 ^,3 76 ._107 _f,2^ E72_A73_i7 72 
2y 'a;m + 27 ~'~ 27 27 9 ~'~ 54 27 9 

2 1 2 A 

(4.6) 

_ 16681 _ 38 8 _HZ7 H72_i73 10921 8 
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^72 7 _^7 73 J_72 72 _J_,3 , /_ 12733 29_ J_ 4 

+ 24 'a:/^ 'a;m -j^g 'i;/:* l-xm + ^2 ^i" Ig + '^M 7775 324 45^ 

^ ^ _ 1 ^ ^ 1^ 72 73 ^74 16^ , 397, 

+ g CS 2 ^3 Ixm + 3ggg Lxm + ^32 ^zm + 215 'xm + 3g ^xm + 27 ^3 Ix^ 324 

2 2 1 2 43 2 1 257 7 2 1 3 

+ 9 Cs ^s^i + Yog ^a:/" ~ 72 ^^1^ 3 ^^'^ ^^"^ ~ 108 ^^"^ 12 ^^'^ ~ 9 

1 2 ,2 _1,3 , A 2f_^^ Za7 ^7 

+ g g Lx^ Lxmj + riii^ ^^^g + 162 90"" ^ 9 ^ 9 ^ 486 """^ 

25 ,2 ^ ;3 J_ 74 _ ^ ^ 7 _ i /- 72 1 ,2 5 , ,2 

+ iQg ''xm + g4 ''xm + 72 'xm 27 '^7* 1296 "'"9 216 18 

~ 108 ^^'^ ^^"^ ~'~ 24 ^^'^ ^^"^ ~ 18 ^^'^ ^^"^ 12 ^^'^ ~ 18 ^^'^ ^^"^^ (4-^) 

_ 205123 593 5 , 15017 , 325 , 47 , 1 4 

^3 " 1049760 2430 ^' 54 "™ + 23328 ^ 3888 5832 ^™ + 324 

^ /- _ 72 J_ 73 _ ^ 7 73 25 ,2 ,2 1 ,3 , 

+ g^C3<x^ + g^4g/xM -^944 + 36 ^o^M 486 '^'^'^"^+ 648 ^'^ ^"* 36 ^a.^'^'" 

109 , 2 _1M^7 ^^,2 49255 167 ^ 7 4 

"1944 ^^"^ " 1944 '""^ + 1944 + V 104976 " 162 + UlE"" 

— r 1^ 7 3149 ,2 175 ,3 7 ,4 ^ , >■ 29663 

+ gj C3 /a.m + 3ggg /xm + jjgg^ ^xm + 5832 ^ 972 ^ 81 ^ 34992 ^'^ 

_^72 , J_73 _ J_ 7 73 , ^ 72 72 _ J_ 73 7 --^7 72 
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1458 972 



2 / 9157 53 ^ 4 4 59 ^ , 5905 , 709 . 7 ,3 



V 209952 243 ^ 1215 162 23328 5832 "^"^ 243 



xm 



1 4 2189 227 2 ^,3 _A/ /3 25 o ,2 

+ 'xm 27 34992 '''' 1944 ^ 36 243 '''' "^"^ ^ 648 ""^^ 

1^ 1 _ / 72 _ 2^2^ 7 / , 72 / ^ (A q^ 

2g 'x/i ''xm 224 'xm 5832 1944 '"^ / \^-^) 

35389447 227 , 895 ^ , 2812987 , 125203 . 4603 

^ ~ 44789760 103680 1728 2985984 248832 ^'"^ 62208 "^'"^ 
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41472 ^'^ ^ 64 '^'^ ^ 165888 '''' 497664 ^ 288 5184 ^''^ ^"^ 

229 ,2 ,2 _ J_ ,3 , _ 8501 2 _ 80177 2279 2 



gg-^2 ^™ 48 '''' 41472 124416 20736 '''' 

^ / 19473727 27725 43 ^4 271 ^ 3637091 226123 2 

' V 35831808 ~ 20736 '"^ ^ 6480^ 1728 ''^ ^™ ^ 2985984 ^"^ 497664 
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_| :_ 73 / /3 _| /2 /2 73 ; ; r2 

288 ^'" 10368 ^^ ^"^ 6912 48 'x/. W 4^472 '^'"^ 



'x/i ^xm ~l~ on^yoc 'xu 'xm I (^'^O) 



35081 2279 o 

"124416 "'^ ""^^ 20736 

with l^f, = log{-p'^ / fi^) and /^.^ = logi-p"^ /m{fif). 

In fig. Q-^ the exact result obtained from numerical integration is compared against the 
asymptotic behaviours of different expansion depths for Minkowskian momenta with m = 
11. In both regimes, we find that already a moderate number of terms reproduces the exact 
curve very well even close to threshold. 



5. Conclusions 

In this work, the low- and high-energy expansions together with an exact numerical solu- 
tion of the double fermionic contribution to the heavy quark vector current correlator in 
four-loop approximation were obtained. We have shown that the method of differential 
equations provides an excellent tool to compute this quantity with high precision in the 
whole momentum region. Consequently, the completion of the IBP reduction represents 
the only remaining task to obtain the full vacuum polarization at this order. 
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Figure 4: Comparison between expansions and numerics of the real part of CpT^^fi' (z) for each 
coefficient of n; separately. Below threshold (represented by the vertical line in z = 1), the dashed 
(solid) curves correspond to low-energy expansions including the first 5 (30) terms. Above threshold 
the dashed, dash-dotted and solid lines denote the high-energy expansions including the first 5, 10 
and 15 terms, respectively. 
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